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ABSTRACT 

We combine imaging from the Hubble Space Telescope Advanced Camera for Surveys, as part of the 
GEMS (Galaxy Evolution from Morphologies and SEDs) survey, with redshifts and rest-frame quantities 
from COMBO- 17 to study the evolution of morphologically early- type galaxies with red colors since 
^ . z = 1. From 0.5° x 0.5° imaging, we draw a large sample of 728 galaxies with centrally-concentrated 

' radial profiles (i.e. n > 2.5 from Sersic fits) and rest-frame {U — V) colors on the red sequence. We 

, explore how the correlations of rest-frame F-band luminosity and of stellar mass with intrinsic half- 

■ light size change over the last half of cosmic time. By appropriate comparison with the well-defined 

local relations from the Sloan Digital Sky Survey, we find that the luminosity-size and stellar mass-size 
relations evolve in a manner that is consistent with the passive aging of ancient stellar populations. By 
' itself, this result is consistent with a completely passive evolution of the red early-type galaxy population. 

If instead, as demonstrated by a number of recent surveys, the early-type galaxy population builds up 
. in mass by roughly a factor of two since z ~ 1, our results imply that new additions to the early- type 

JZLh' galaxy population follow similar luminosity-size and stellar mass-size correlations, compared to the 

older subset of early- type galaxies. Adding early- type galaxies to the red sequence through the fading of 
^ ■ previously prominent disks appears to be consistent with the data. Through comparison with models, 

^ ' the role of dissipationless merging is limited to < 1 major merger on average since z = 1 for the most 

, massive galaxies. Predictions from models of gas-rich mergers are not yet mature enough to allow a 

detailed comparison to our observations. We find tentative evidence that the amount of luminosity 
evolution depends on galaxy stellar mass, such that the least massive galaxies show stronger luminosity 
evolution compared to more massive early types. This could reflect a different origin of low-mass early- 
' type galaxies and/or younger stellar populations; the present data is insufficient to discriminate between 

5^ I these possibilities. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters (luminosities, stellar masses, 
radii) — galaxies: general — surveys 



1. INTRODUCTION 

The formation and evolution of massive, early- type galax- 
ies constitutes a long-standing and crucial problem in cos- 
mology. In all hierarchical models, e.g. a A-dominated 
Cold Dark Matter (ACDM) cosmology, the massive early- 
type galaxies seen now are expected to have formed through 
mergers of smaller galaxies over time (White & Frenk 1991; 
Cole et al. 2000). Yet, within these models it is difficult to 
predict robustly at what point during this assembly most 
stars formed. In order to constrain the star formation (SF) 
and assembly histories of the early-type galaxy population, 
it is necessary to explore the evolution of their number 
density and various scaling relations, such as the funda- 
mental plane (van Dokkum & Franx 1996; Kelson et al. 
2000; van Dokkum et al. 2001; Treu et al. 2002; Wuyts 
et al. 2004; van der Wei et al. 2004), the luminosity-size 



relation (Barrientos et al. 1996; Pahre et al. 1996; Schade 
et al. 1997; Barger et al. 1998; Schade et al. 1999; Ziegler 
et al. 1999), and the stellar mass-size relation (Trujillo 
et al. 2004b). To this end, we use the GEMS (Galaxy 
Evolution from Morphology and SEDs; Rix et al. 2004) 
survey in conjunction with COMBO-17 (Classifying Ob- 
jects by Medium-Band Observations in 17 Filters; Wolf 
et al. 2003, 2004) to construct the largest sample to date 
of distant (0.2 < z < 1.0) early-type galaxies with Hubble 
Space Telescope (HST) imaging to explore the evolution of 
luminosity and stellar mass as functions of galaxy size. 

Most initial studies of z < 1 early-type galaxies focused 
on galaxies in some of the richest clusters at each epoch 
(e.g., van Dokkum & Franx 1996; van Dokkum et al. 1998; 
Kelson et al. 2000), and found that the mass-to- light ratios 
(M/L) of massive galaxies in this environment changed as 
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expected for passive fading of a population formed exclu- 
sively at very early epochs. In recent years, greatly ex- 
panded studies of "red" galaxies in random cosmological 
volumes have been carried out, extending to redshifts of at 
least z ~ 2 (Chen et al. 2003; Rusin ct al. 2003; Bell ct al. 
2004b; Drory et al. 2004; Fontana et al. 2004; McCarthy 
et al. 2004). From colors and spectra, these studies in- 
ferred that massive galaxies existed at all observed epochs 
and that red galaxies at any given epoch contain almost 
exclusively old (relative to the observational epoch) stel- 
lar populations. Yet, the age distribution of the stars in 
these galaxies does not by itself speak to their dynamical 
assembly. 

While the optical colors of massive red early-type galax- 
ies would indeed be consistent with a population that has 
evolved only through fading since z-^ 2, other data in- 
dicate that more must have happened. In particular, the 
co-moving total stellar mass density in red-sequence galax- 
ies is lower at earlier epochs, amounting to a factor of 
~ 2 buildup since z ~ 1 (Chen et al. 2003; Bell et al. 
2004b; Cross et al. 2004). Several factors are expected 
to contribute to this evolution: (1) in galaxies with wan- 
ing and disappearing SF, the disk will fade and the bulge 
gain prominence in comparison; (2) spheroidal components 
may newly form from gas-rich disk mergers or from in- 
ternal disk instabilities; or (3), galaxies within the red 
sequence could undergo dissipationless merging, changing 
their masses and structures while; Iciaving their stellar pop- 
ulations unchanged. All processes are expected to con- 
tinue to the present epoch in a hierarchical universe (e.g., 
Cole et al. 2000; Steinmetz & Navarro 2002; Khochfar & 
Burkert 2003). The relative importance of the effects is 
expected to be a strong function of galaxy mass: gas-rich 
mergers, disk instabilities, and disk fading are expected to 
contribute strongly at low masses, whereas dissipationless 
mergers may dominate for high-mass galaxies (Khochfar 
& Burkert 2003). 

The aging of a stellar population not only implies red- 
der colors but also dimming. We exploit this to apply an 
independent test for the population evolution of massive 
galaxies. This test uses the morphological and structural 
information for a large sample of early-type, red galaxies 
to z 1, which we draw from the COMBO- 17 redshift 
and spectral energy distribution (SED) survey (Wolf et al. 
2003, 2004) and the GEMS ffST imaging survey (Rix et al. 
2004). By constructing the luminosity size (Ly — rso) and 
stellar mass-size (M* — r^o) relations for subsamples of 
morphologically selected early-type, red galaxies at differ- 
ent redshifts, one can test the hypothesis of passive evo- 
lution and provide constraints for other evolutionary sce- 
narios. In practice, the concentration of the light profile, 
specified by the Sersic index n > 2.5 is used as a quanti- 
tative proxy for morphology. 

For the most strict version of passive evolution (no new 
stars; no merging) the prediction is clear: the — r^Q 
relation will remain unchanged ^ and the Ly — r^o rela- 
tion should change in the sense that galaxies of a given 
size would have been brighter by about 1 magnitude in V- 

^ This is true neglecting the possible change in effective radius that 
might occur in a galaxy with a substantive age gradient. Stellar mass 
loss in an aging stellar population would lead also to slow changes 
in the M* — r^Q relation. 



band at z = 1, if the observed color evolution (Bell et al. 
2004b) is a guide to the stellar M/L evolution. The pre- 
dictions for the evolution of these relations in the presence 
of merging are maturing rapidly. Gas-free dissipationless 
mergers (sec e.g., Naab & Burkert 2003; Dantas et al. 2003; 
Nipoti et al. 2003) will slowly move galaxies away from the 
Ly — r5o and M,i, — r^o relations (Nipoti et al. 2003); we 
discuss this issue in detail later. Existing constraints from 
the color-magnitude relation (CMR) or fundamental plane 
(FP) limit the amount of dissipationless merging that can 
occur to a factor of a few in terms of mass growth in this 
merging mode, if roughly equal-mass mergers dominate 
the mass growth (Bower et al. 1998; Dantas ct al. 2003; 
Nipoti et al. 2003; Gonzalez-Garcfa & van Albada 2003; 
Evstigneeva et al. 2004). Predictions for gas-rich mergers 
are less mature, and it is not clear what to expect for the 
evolution of the Ly — r^o and M,^ — r^o relations in this 
case. 

As shown persuasively by Simard et al. (1999), any such 
test for passive or other physical evolution requires care- 
ful accounting for redshift-dependent selection effects. We 
will address this issue extensively, as the systematics of 
the sample selection become increasingly important over 
diminishing statistical errors for larger samples. Also, all 
existing surveys that are deep enough to reach z ^ 1 
cover such a small area that their volumes at low redshifts 
(z < 0.2) become too small. Therefore, we will combine 
the C0MB0-17/GEMS data at z > 0.2 with the present- 
day Ly — r5o and stellar M* — rso relations for early-type 
galaxies (Shen et al. 2003) from the Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000). Both for Shen et al. (2003) 
and the 0.2 < z < 1.0 sample presented here, stellar mass- 
to-light ratios are estimated from the spectra and SEDs. 

In this paper, we bring to bear the largest sample of mor- 
phologically and color-selected early-type galaxies with HST 
imaging for a comprehensive analysis of the evolution of 
the Ly — rso and M,^ — rso relations. Measuring the evo- 
lution of relationships between these fundamental galaxy 
properties over < z < 1 can place further constraints 
on galaxy formation and evolution models. Specifically, 
we have two goals: (1) to quantify the expected luminos- 
ity evolution of galaxies of a given size due to the simple 
aging of their ancient stars with the best accounting of se- 
lection effects to date; and (2) to test whether stellar mass 
evolution occurs at different fixed galaxy sizes, which may 
provide constraints on the merging history of these systems 
of old stars. The outline of this paper follows. We describe 
our early- type galaxy data in §2, including relevant discus- 
sions of the GEMS imaging survey, the rest-frame quanti- 
ties from COMBO-17, the early-type sample selection, and 
our galaxy size measurements. In §3 we explore the evolu- 
tion of the luminosity and stellar mass scaling relations of 
early-type galaxies with z < 1. We discuss our results in §4 
and give our conclusions in §5. For rest-frame luminosi- 
ties, stellar masses, and physical sizes we assume a flat, 
A-dominated cosmology with = 0.3, Qa = 0.7, and 
Ho = 100/ikms~^ Mpc~^. When estimating expected pas- 
sive stellar population evolution, we adopt a flat universe 
with the WMAP parameters {Hq = 72kms"i Mpc^^ and 
Qui = 0.27, Spergel et al. 2003), which give a Hubble time 
of 13.5 Gyr, to account for stellar M/L and color evolution 
in the most realistic way possible at this time. 
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2. EARLY- TYPE GALAXY DATA 

We select a well-defined sample of early-type galaxies 
with i/ST imaging from GEMS, and high-precision photo- 
metric redshifts and rest-frame luminosities from COMBO- 
17. In this section we discuss briefly the high-resolution 
imaging and the aspects of the COMBO- 17 data relevant 
to this study. Furthermore, we outline the early-type sam- 
ple selection and its completeness, and wc describe the 
galaxy size measurements from the HST data. 

2.1. High-resolution Imaging 

The GEMS survey has obtained a large, two passband 
(F606W and F850LP) Advanced Camera for Surveys (ACS) 
image mosaic over an area of 28' x 28', encompassing the 
Extended Chandra Deep Field South (E-CDFS). This is 
the largest contiguous color map of the cosmos obtained 
with HST and consists of a grid of 78 mostly overlap- 
ping images taken with the ACS wide field camera dur- 
ing November 2002. The full details of the imaging grid, 
observations, data reduction, calibration, and data quality 
assessment will be presented in Caldwell et al. (in prepara- 
tion, C05). An overview of the GEMS survey experimental 
design is given in Rix et al. (2004). The foremost goal of 
GEMS is to quantify the internal structural evolution of 
galaxies using statistically significant samples of galaxies 
at known redshifts. To this end we selected the E-CDFS 
field with an existing redshift data set from COMBO- 
17. Furthermore, a fraction of GEMS overlaps with the 
GOODS (Giavahsco et al. 2004), Ultra Deep Field (UDF), 
and UDF parallels (Bouwens et al. 2004) programs, thus 
providing a wealth of deeper, multi-wavelength data. 

Hereafter, wc concentrate our analysis on measuring 
galaxy sizes from the longest rest-frame wavelength possi- 
ble at which galaxy profiles appear most uniform and the 
contamination from star-forming regions in disk galaxies 
is minimized. Therefore, we measure sizes by fitting sim- 
ple models to the two-dimensional galaxy imaging in the 
F850LP filter. 

For each pointing in the GEMS mosaic we have indi- 
vidual F850LP images, which are multi-drizzle combina- 
tions of three dithered exposures (720-762 seconds each; 
see C05 for details). Briefiy, the exposures were processed 
individually prior to combining (i.e. bias and dark cur- 
rent correction, flat fielding). The final frames are free 
of most cosmic rays and remapped onto a fine O'.'03/pixel 
scale. The full suite of 78 F850LP images are quite flat 
with background differences between images considerably 
less than 0.2 ADU, and the local galaxy background mea- 
surements across each frame have a typical RMS of 0.2 
ADU. In addition, wc have a variance image (necessary for 
source detection) with the same scale for each frame. The 
images are flux calibrated for photometric uniformity over 
the large mosaic of ACS frames, which assures unbiased 
galaxy size measurements. We calibrate the astrometry 
of each image to the ground-based epoch J2000.0 system 
of COMBO-17, which allows precise image/redshift cross 
correlation. The reduced frames have angular resolution of 
0'.'077 FWHM, corresponding to a physical scale of 700 pc 
at z ~ 0.75, comparable to the resolution of Coma galaxy 
cluster observations with 1" seeing. 

The final GEMS source catalog contains 41,681 unique 
objects detected in the F850LP imaging. For source de- 



tection we employ SExtractor (Bertin & Arnouts 1996) 
in a two-step strategy discussed in detail in Rix et al. 
(2004). Briefly, flrst we use a conservative ("cold") de- 
tection and deblending conflguration that avoids spurious 
deblending of large galaxies with strong features, yet is in- 
complete for faint, low surface brightness objects found in 
COMBO-17 . Next we use a "hot" configuration to detect 
99% of known COMBO-17 sources down to a total appar- 
ent magnitude of i?tot = 24 (Vega), the limit at which 
the redshift performance drops dramatically. The exact 
configuration parameters are given in Rix et al. (2004). 
In each case, separate catalogs for each GEMS pointing 
arc constructed from the combined F850LP image and a 
weight map (oc var~^). The use of weight maps reduces 
the number of spurious detections in low signal-to-noise 
(S/N) areas of each image (e.g., near image edges). 

For each field we combine all "cold" sources plus the 
subset of "hot" sources residing outside of the isophotal 
area of any cold detections, and then we remove dupli- 
cate sources to produce the final catalog with the following 
properties: (1) contains 99% of all objects in the i?tot < 24 
mag COMBO-17 sample; (2) avoids spurious deblending 
of large, bright galaxies exhibiting strong substructure in 
HST images from spiral features, etc.; and (3) provides a 
homogeneous, flux and surface brightness limited catalog 
of all sources in the F850LP ACS mosaic, regardless of 
COMBO-17 or other external information. Through de- 
tailed simulations (Haufiler ct al., in preparation, hereafter 
H05) we find that galaxy detection in GEMS is complete 
to effective surface brightnesses in the F850LP-band of 
^ 25 mag arcscc^^ (within the half-light radius) for 
r^/* spheroids and /iso ^ 24 mag arcsec^^ for exponen- 
tial disks. The GEMS SExtractor source catalog will be 
published and described in complete detail in COS. 

2.2. Rest-frame Quantities from COMBO-17 

COMBO-17 is a deep photometric survey comprised of 
five (three completed) disjoint fields, each ^ 0.25 square 
degree in size, with flux measurements in 5 broad and 12 
medium passbands (between 3500 and 9000 angstroms). 
This survey provides data for ~ 25, 000 galaxies with iitot < 
24 mag in the range 0.2 < z < 1.2. The COMBO-17 data, 
associated uncertainties, sample selection, and complete- 
ness are discussed in detail in Wolf et al. (2004). 

The deep 17-passband data combined with stellar, AGN, 
and non-evolving galaxy template spectra allow nearly all 
(98%) objects to be assigned redshifts and SED classifica- 
tions. The redshift accuracies depend primarily on i?-band 
aperture magnitude, with 5z/{l-\- z) ~ 0.02 for iiap < 22 
mag. When going fainter the errors reach 5z/{\-\-z) ^ 0.05 
at i?ap = 24 mag (Wolf et al. 2004). Reliable redshift es- 
timates allow the construction of rest-frame luminosities 
and colors. Depending on the redshift in question, rest- 
frame luminosities are determined either by interpolation 
or mild extrapolation. For this study we use the rest-frame 
V"-band absolute magnitudes for the luminosity parameter, 
which are extrapolated for galaxies with z > 0.7. Never- 
theless, the leverage of the COMBO-17 filter set continues 
to provide accurate F-band magnitudes to 2; = 1, with typ- 
ical uncertainties of 10% (z < 0.7) and 15%. (z > 0.7). 
We note that the quoted redshift accuracies correspond to 
< 15% fractional distance errors for z > 0.2. 
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We cross correlate the positions of sources from the 
GEMS F850LP imaging with source coordinates from the 
COMBO-17 E-CDFS catalog^. Wc find 6152 galaxies with 
redshifts 0.2 < z < 1.0 and i?ap < 24 mag. This is our 
overall galaxy sample. Wc note that the COMBO-17 mag- 
nitude limit corresponds to the reliability of estimating 
redshifts, which is brighter than the limit to detect galaxies 
with high completeness both in COMBO-17 and GEMS. 
For red-selected galaxies, which are the primary concern 
of this paper, the completeness limit owing to redshift re- 
liability depends somewhat on redshift: for 0.2 < z < 0.6 
the mean 90% limit occurs at R^p ~ 23.5, while at larger 
redshifts (0.6 < z < 1.0) the sample is fully complete down 
to i?ap = 24 mag. At redshifts above z = 1, COMBO-17 
remains complete to the magnitude limit, yet this limit 
begins restricting the sample to only the very luminous 
(My < —21 + 5 logiQ /i at z ~ 1) red galaxies. We there- 
fore adopt a strict z = 1 cut for our analysis. 

We use stellar mass M^, estimates based on direct mod- 
eling of the COMBO-17 SEDs (Borch et al., in prepara- 
tion). The 17-passband COMBO-17 SEDs were compared 
with a library of galaxy template spectra obtained from 
the Pegase code (sec Fioc & Rocca-Volmcrangc 1997, for 
a description of an earlier version of the code). From each 
galaxy's SED we estimate directly its stellar M/L ratio 
(M*/L), and hence M*; M* uncertainties include contri- 
butions from redshift uncertainty. We use a Kroupa et al. 
(1993) stellar initial mass function (IMF), which produces 
Mi,/L values comparable to those from a Kroupa (2001) 
IMF, and 0.3 dex lower than from a Salpeter (1955) 
IMF. For rcd-sclcctcd galaxies, these stellar masses agree 
to within < 40% (0.11 dex), and show no average offset 
when compared to stellar masses based on the simple re- 
lation between {B — V) color and M*/L from Bell et al. 
(2003), assuming the same IMF. The stellar mass esti- 
mates suffer from a number of random and systematic un- 
certainties, such as uncertainties in galaxy age, bursts of 
SF in the last 1-2 Gyr, and variations in the relationship 
between dust reddening and extinction (see e.g.. Bell & de 
Jong 2001; Bell et al. 2003, for a more in-depth discussion 
of the relevant sources of uncertainties). Generally these 
uncertainties amount to ~ 0.3 dex; for the red-sequence 
galaxies studied here, it is likely that the stellar mass esti- 
mates are more accurate than this, as relatively old stellar 
populations are somewhat easier to model robustly with 
simplistic SFH prescriptions. 

2.3. Early-type Galaxy Selection 

In very generic terms, the visual appearances of galaxies 
fall into two broad types: (1) early- type refers to a galaxy 
with spheroid or bulge-dominated (centrally-concentrated) 
morphology; and (2) late- type corresponds to a disk-dominated 
or irregular system. In the local universe, several classes 
from the familiar Hubble (1926) sequence are considered 
early- type in morphology - ellipticals (E), lenticulars (SO), 
and spirals with dominant bulges and small disks with 
tightly wound spiral arms (Sa). In addition to physical 
appearance, the optical colors of galaxies are broadly dif- 
ferent - early types have typically red colors suggestive of 
a dominant population of old stars with little or no current 

^ Available at littp;//www.mpia- 

hd.mpg.de/COMBO/comboJndex.html. 



SF; late types tend to be bluer in color as a result of on- 
going SF and some fraction of young stars. This bimodal 
distribution of galaxies in color space is observed both lo- 
cally (Strateva et al. 2001; Hogg et al. 2002; Blanton et al. 
2003b; Baldry et al. 2004) and out to z - 1 (Bell et al. 
2004b; Weiner et al. 2004). Furthermore, red galaxies at 
z ^ 0.7 contain the same predominance (^ 85% by op- 
tical luminosity density) of early- type galaxies (E/SO/Sa) 
as locally (Bell et al. 2004a), with a small contamination 
by edge-on disk galaxies with colors reddened by internal 
extinction. For this study we are interested in early-type 
galaxies with red colors, thus we select our sample using 
both color and morphology. 

2.3.1. Red Sequence 

Red-color selection permits an objective and empirical 
first cut for defining early-type galaxies over the last half of 
cosmic history. We use the empirical fit to the CMR evo- 
lution (including local comparison points) from Bell et al. 
(2004b)^ to select red-sequence galaxies at a given redshift 
z with rest-frame colors redder than 

{U-V) = 1.05 - 0.31z - 0.08(My - 5 logip h + 20). (1) 

This color cut is 0.25 mag redder than the fit to the CMR 
evolution, and similar in philosophy to the Butcher & 
Oemler (1984) criteria, except that wc select red rather 
than blue galaxies. The color evolution of red-sequence 
galaxies selected with this definition is in good agreement 
with the expected evolution of a single-age stellar popu- 
lation anchored to z = 0. There are 1166 red-sequence 
objects with 0.2 < z < 1.0 and GEMS imaging. 

To facilitate meaningful profile fitting, we visually in- 
spect the F850LP images of the red-sequence sample. COMBO- 
17 is expected to misclassify a small number of red M-class 
stars as modcrate-rcdshift red-sequence galaxies (Wolf et al. 
2004). With the superior resolution of GEMS, we can di- 
rectly detect 69 misclassified stars, which are shown in 
Figure 5 for completeness, and are omitted from the anal- 
ysis hereafter. We note that these stars are also separated 
cleanly from galaxies using a SExtractor CLASS_STAR> 
0.85 cut. Additionally, we remove 50 red-sequence galaxies 
with poor S/N found near GEMS image edges or within 
inter-chip gaps. This leaves 1047 red-sequence galaxies 
with useful GEMS imaging. Last, we remove 84 galax- 
ies with morphologies where reliable fits are not possible 
(14 with prominent dust lanes, 29 irregulars, and 41 pe- 
culiar/interacting). We are left with a final sample of 963 
morphologically "normal" red galaxies with the following 
breakdown of coarse types: 699 E/SO, 162 Sa, and 102 
late- type spirals (Sb-Sd). Owing to the subjective nature 
of visual classifications, we divide the sample morphologi- 
cally using a quantitative method in the next section. 

2.3.2. Quantitative Morphology 

We refine our sample further by determining which red- 
sequence galaxies have early-type morphology using a quan- 
tifiable and repeatable method. The simple Sersic (1968) 

^ This cut is slightly bluer than that given by Bell et al. (2004b). Ow- 
ing partially to uncertainties in photometric calibration and galactic 
foreground extinction, the red sequence has slightly different colors 
in the three COMBO-17 fields, and the cut we use here is more 
appropriate for the E-CDFS. 
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model, which describes the radial surface brightness pro- 
file of a galaxy by I](r) oc exp (— r-^/"), is the most widely- 
used parametric function. The Sersic index n describes 
how centrally-concentrated a galaxy appears, with n = 4 
describing the familiar r^^^ profile of ellipticals (dc Vau- 
couleurs 1948), and n = 1 representing the exponential 
profile commonly seen in spiral disks. Studies with the 
SDSS have adopted empirical cutoffs around n = 2.5 to 
separate morphologically early and late types (Blanton 
ct al. 2003b; Shcn ct aL 2003; Hogg ct al. 2004). In ad- 
dition, using a sample of nearly 1500 GEMS galaxies in 
the thin redshift slice 0.65 < z < 0.75, Bell et al. (2004a) 
showed that an n = 2.5 cut was reliable at distinguishing 
between visually-classified early and late types. To remain 
consistent, we adopt n > 2.5 to denote early types in this 
study. 

We note that the n = 2.5 cut is a crude concentra- 
tion cut that takes no account of the real diversity of 
galaxy forms seen in nature. In particular, this cut does 
not account explicitly for dynamically-distinct bulge and 
disk components, and we reserve a detailed fitting of the 
photo-bulge and photo-disk components to future papers. 
Nevertheless, it is unclear whether bulge/disk decomposi- 
tions will be any more robust or physical than the simple 
n = 2.5 cut. We note two additional sources of concern 
with the n > 2.5 selection. First, (1 + 2;)^ dimming of disks 
may cause a systematic bias such that large disks at high 
redshift may bo too faint to detect and, thus, appear as 
bulge-dominated galaxies (larger n values). Second, a cor- 
relation exists between n and absolute B-band magnitude 
for visually-selected ellipticals (e.g., Trujillo et al. 2004a). 
Therefore, a strict n = 2.5 cut is effectively a luminosity 
cut of roughly My = —17+5 log^Q h. This magnitude Hmit 
is below our redshift reliability cutoff at all redshifts ex- 
cept for the z < 0.3 interval, at which the drop in numbers 
of galaxies fainter than —17.5 + 5 logj^Q h may be explained 
by this effect. Ultimately, we accept these caveats with the 
understanding that we have a repeatable and empirically- 
motivated early-type galaxy sample. 

In Appendix A, we describe our method for fitting a 
Sersic model to each galaxy's two-dimensional luminosity 
profile from the F850LP imaging. Briefly, we use the fit- 
ting code gim2d (Simard et al. 2002) to fit the sample of 
963 red-sequence galaxies and find that 96% (928) are fit 
successfully. The majority of fitting failures are the result 
of GIm2d reaching the Sersic parameter limit of n = 8. We 
estimate uncertainties of 5n/n ~ 0.25. Our final sample 
contains 728 red galaxies out to z = 1 with n > 2.5. In 
general, using an n = 2.5 cut does a good job separating 
red-selected galaxies into morphologically early and late 
types. We find only 3% of the red n > 2.5 sample have 
late-type (Sb-Sd) visual morphologies, while 14% of the 
visually early types (E-Sa) have n < 2.5; this is typical 
of what was found in the comparison of by-eye and Sersic 
classifications a.t z ^ 0.7 by Bell et al. (2004a). We show 
in Figures 1-3 that representative examples of the red 
GEMS galaxies with n > 2.5 are visually early- type in ap- 
pearance. For contrast, in Figure 4 we display a handful 
of the red-sequence galaxies with n < 2.5 that correspond 
to edge-on disks. 

2.4. Effective Size Measurements 



2.4.1. Galaxy Sizes 

The main focus of this work is the measurement and 

analysis of galaxy sizes as a function of redshift. GIm2d 
calculates the half-light semi-major axis length a^o by in- 
tegrating the best-fit Sersic model flux to infinity. For 
meaningful comparison with SDSS results we adopt a cir- 
cular half-light radius, or geometric mean, given by r^Q = 
tJsovl — e, where e is the best-fit model ellipticity. We 
have applied a robust method for determining the system- 
atic and random uncertainties of each model parameter by 
measiiring their random offset and variance with respect 
to known values using a large sample of artificial galaxies 
analyzed in the same manner as the GEMS data (H05). 
This method has proven successful in other galaxy fitting 
work (e.g., Tran et al. 2003; MacArthur et al. 2003). We 
find that rso has typical random imcertainties of ^ 35% 
with no systematic error. In addition, we find that GIm2d 
r^o measurements arc reliable for > 90% of galax;ies with 
F850LP-band surface brightness of /iso < 24 mag asec~^. 
As seen in Figure 5, our early- type sample is well within 
this limit. 

Sky level uncertainties are an important source of sys- 
tematic error in size measurements when applying profile 
fitting techniques (e.g., de Jong 1996). We estimate the 
error in our fitting-derived sizes due to the micertainty in 
our local background measurements ((Tbkg ^0.2 ADU) by 
repeating the fits with the sky level held to constant val- 
ues of ilcTbkg- We find that our size measurements have 
average uncertainties due to the sky of ^ 10%. 

Finally, our galaxy sizes are based on r^Q in a fixed, ob- 
served passband (F850LP), which corresponds roughly to 
the SDSS 2;-band. Therefore, galaxy sizes at different red- 
shifts are measured at different rest-frame wavelengths; 
e.g., the central wavelength of the F850LP filter corre- 
sponds to rest-frame 7560 Aat z = 0.25 and 4846 Aat z = 
0.95. For a consistent comparison with the local scaling 
relations from Shen et al. (2003), we correct each rso mea- 
surement to rest-frame r-band. Early-type galaxies are 
known to have internal radial color gradients (e.g., Franx 
et al. 1989; Peletier et al. 1990; Bernardi et al. 2003), with 
redder colors towards their centers. Such gradients should 
result in a wavelength-dependent galaxy size. We esti- 
mate a correction to shift z-band sizes to rest-frame r-band 
based on typical color gradients of A{B — R)/A\ogr = 
-0.09 and A{U - i?)/Alogr = -0.22 from Peletier et al. 
(1990). We assume an r^^^ surface brightness profile and 
find that this gradient will produce measured sizes in the 
[/-band to be about 6.5% bigger than in the r-band; i.e. 
Ar/AlogA ~ —0.25. Between z = 1 and z = the same 
observed passband corresponds to a factor of 2 change in 
wavelength, or a 7.5% size correction per unit redshift. 
The correction to rest-frame r-band (^-band observed at 
z ~ 0.45) is then given by 
Ar 

= -0.075(z- 0.466). (2) 

^obs 

Owing to the mild color gradients this correction is quite 
small, amounting to only 4.0% (1.5%) decrease (increase) 
in observed size for z = 1 (z = 0.25), and as such it does 
not significantly affect our conclusions. We note that ellip- 
ticals with blue cores (i.e. inverted color gradients) have 
been observed by Menanteau et al. (2001) in the Hubble 
Deep Field (HDF), yet these galaxies are globally bluer 
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than our red-sequence cut. The assumptions we make re- 
garding color gradients are based on red galaxies; there- 
fore, it is fair to apply this correction to our red-selected 
sample. 

2.4.2. Magnitude- Size Distributions 

Before proceeding with our analysis of scaling relations 

using rcst-framc quantities, we plot in Figure 5 the magnitude- 
size distributions in the observed F850LP-band frame for 
all 928 red-sequence galaxies (split into eight bins of width 
Az = 0.1 spanning redshifts 0.2 < ^; < 1.0) to illustrate 
several aspects of our sample. We calculate the apparent 
magnitude mgso of each galaxy from the total intensity of 
the best- fit model (see Appendix A). We see that galaxy 
apparent magnitudes and angular sizes are correlated out 
to z = 1, such that bigger galaxies are brighter. We note 
that the effects of (1 + z)^ cosmological dimming appear 
as a shift of the magnitudo-sizc correlation towards lower 
surface brightnesses at increasing redshifts. Specifically, 
between z = 0.25 and z = 0.75 the mean relation has faded 
by the expected 1.5 mag. Furthermore, nearly all galax- 
ies have rgo < I'.'O, clearly illustrating the need for HST 
resolution to obtain accurate size measurements for dis- 
tant galaxies. All red-selected galaxies are well-resolved, 
and have effective surface brightnesses > 1.0 mag brighter 
than the GIm2d reliability limit and > 2.0 mag brighter 
than the 85% detection completeness limit. Therefore, our 
early-type sample is limited only by the i?,ap = 24 mag cut- 
off (corresponding to mgso ^ 22) for reliable photometric 
redshifts inherent to the COMBO-17 survey. 

In addition, we show where n < 2.5 (open squares in Fig- 
ure 5) red-selected galaxies arc found in the magnitude-size 
plane. Most red n < 2.5 galaxies tend to have lower surface 
brightnesses relative to the spheroid- dominated n > 2.5 
galaxies (solid symbols), and are for the most part highly- 
inclined disk galaxies. These results are in qualitative 
agreement with the different magnitude-size correlations 
for reddened edge-on disks and early-type galaxies in the 
local universe (Blanton et al. 2003b). 

3. SCALING RELATIONS OF EARLY- TYPE GALAXIES WITH 

0.2 < z< 1.0 
3.1. Luminosity-Size Relation 

We now present the Ly — r^o relations for early-type 
galaxies with 0.2 < z < 1.0 and we will compare these to 
the present-day distribution. In Figure 6, we divide our 
sample into eight redshift bins of Az = 0.1 and for each we 
plot the absolute F-band magnitude against the rest-frame 
r-band size in physical units (/i^^kpc). The bulk of the 
early- type population at each epoch spans roughly a 0.5 to 
5.0 /i~^kpc half-light size distribution, with a handful of 
smaller (< 0.5/i~^kpc) and larger (> 5.0/i^^kpc) galax- 
ies in some redshift slices. Starting at low redshift, the 
average number of galaxies per redshift slice increases as 
we sample larger co-moving volumes, until z ^ 0.75 where 
the counts start falling off as expected in a magnitude- 
limited sample. The second lowest redshift bin is quite 
sparsely-populated as a result of large-scale fluctuations. 

For comparison, in each redshift interval of Figure 6 we 
show the "ridgeline" of the present-day Ly — rso relation 
for early- type galaxies. The ^; ~ ridgeline is given by 
the median absolute F-band magnitude as a function of 



galaxy size using a mock catalog of My — 5 logj^g ^ 
rso values that follow the luminosity function of SDSS- 
selected, red-sequence galaxies from Bell et al. (2004b), 
and the Ly — rso relation for n > 2.5 galaxies in SDSS from 
Shcn et al. (2003). Using a complete sample of n > 2.5 
galaxies'* these authors found a Ly — rso scaling relation 
in the local universe that is well fit by a simple power 
law, and they showed that SDSS early types of a given 
luminosity follow a log-normal size distribution given by 



f{r50,r5o{L),a{L)) 



1 



27r£7(L) 



exp 



ln^(rso/rso(-J')) 



2(72 (L) 



dri 



50 



r-50 



(3) 

with characteristic median size rso(i) and dispersion ct(L). 
The Shen et al. (2003) work is based on half-light radii and 
total absolute magnitudes from Sersic (PSF-corrected) fits 
to the azimuthally-averaged surface brightness profiles in 
the r-band (for radial profile fitting see Stoughton et al. 
2002). These present-day (0.05 < z < 0.15) sizes arc given 
in observed r-band; nevertheless, at a median redshift of 
z ^ 0.1 the sizes are within 1% of r-band rest-frame ac- 
cording to our passband correction calculation (§2.4). Us- 
ing fc-corrections from Blanton et al. (2003a), the SDSS 
galaxy luminosities have been corrected to the r-band rest- 
frame (z = 0), which we translate to rest-frame 1^-band in 
our adopted {h is free parameter) cosmology using a mean 
{V -r) = 0.33 color for E/SOs (Fukugita et al. 1995). 
Therefore, as a function of absolute U-band magnitude 
My = My — 5 logiQ h, the local Ly — rso relation has me- 
dian r-band size (in units of kpc) 

f50(M{,) = 10(-0-26M^-4-93)^ (4) 

with dispersion^ that increases with decreasing luminosity 

as 



CT(My) = 0.27 + 



0.18 



(5) 



Equation (3) does not directly allow one to plot the 
ridgeline of the Ly — rso relation, because it does not 
account for the steeply-declining luminosity function of 
early-type galaxies, which at large size will shift the ridge- 
line towards fainter luminosities. To account for this we 
construct a large grid of galaxy luminosities and sizes that 
sample the Ly— rso relation of 2; ~ early types with a log- 
normal dispersion in the size direction given by Equations 
(3-5). We then populate this grid so that the Afy — 5 log^o h 
distribution matches the local S-band luminosity function 
given in Bell et al. (2004b), which we transform to rest- 
frame F-band assuming a typical {B — V) = 0.9 color 
for E/SO types (Fukugita et al. 1995). The grid spans 
-24 < My - Slogio h < -16 in 800 (0.01 mag) cells and 
< rso < 25 /j-i kpc in 500 (0.05 h''^ kpc) cells, for a fine 
grid of 400,000 total cells with over 16.9 million Ly,r5o 
values representative of the Ly — rso distribution of early- 
type galaxies in the present-day universe. We illustrate 
this mock catalog by plotting its median My — 51og]^o''- 

Wc note that the Shcn ct al. (2003) carly-typc sample is not 
color selected, thus, there will be some contamination from centrally- 
concentrated galaxies with blue colors. Yet, Shen et al. have stated 
that the majority of n > 2.5 SDSS galaxies have red colors of 
(g ~r) > 0.7 (see also Blanton et al. 2003b). 

^ This expression fits the dispersion of the Ly — rso relation for both 
early and late types to within their respective error bars (see Fig. 6 
of Shen et al. 2003). 
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Fig. 1. — Representative examples of small (fixed-size bin I: 0.5 < rso < 1.0 h^^ kpc) early-type {n > 2.5 profiles with red-sequence colors) 
galaxies in the redsliift range 0.2 < z < 0.7. All galaxies shown arc visually early types. In each row we show the F850LP-band galaxy and 
Gim2d fitting residual postage stamp images for three randomly-selected galaxies in each redshift interval. We give the galaxy COMBO-17 
identification number for each galaxy in the upper left. The postage stamp images axe 7h~^ kpc per side. The horizontal white lines represent 
1 arcsecond. 
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Fig. 2. — Representative examples of medium size (bin li: 1.0 < rgo < 2.0 h~^kpc) early-type galaxies with n > 2.5 profiles and red 
colors spanning redshifts 0.2 < z < 0.9. All galaxies shown are visually early types. Three example images with their corresponding fitting 
residuals are shown from each redshift interval as in Figure 1. The postage stamp images are 14h~^ kpc per side, with 1 arcsecond denoted 
by horizontal white lines. 
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Fig. 4. — Examples of late-type disk galaxies (n < 2.5 profiles) with red-sequence colors. We show the F850LP-band galaxy and GIM2d 
fitting residual postage stamp images for 8 edge-on disks spanning rcdshifts 0.25 < z < 0.73. We give the galaxy COMBO-17 identification 
number and redshift for each galaxy. The postage stamp images are I4h~^ kpc per side. The horizontal white lines represent 1 arcsecond. 



(with 16-percentile and 84-percentile limits) as a function 

of r-band size in the left panel of Figure 7. As can be 
seen in Figure 6, the Shen et al. (2003) ridgcline provides 
a reasonable fit to the GEMS Ly — r^o relations out to 
z 0.5. At larger redshifts the relations start to deviate 
from the local power law, with the largest evolution ap- 
pearing for the smallest galaxies. In any galaxy evolution 
scenario we expect that passive luminosity evolution of the 
stellar populations in galaxies must play a role. We use 
the mock catalog to study the evolution of the Ly — r^o 
relation with redshift in the next section. 

3.2. Luminosity Evolution at Fixed Galaxy Sizes 

We now quantify the evolution seen in Figure 6, which 

can be interpreted in terms of passive evolution of an an- 
cient stellar population. At a given fixed galaxy size, we 
measure the luminosity difference relative to the z ~ 
relation as a function of redshift, accounting for selection 
effects. With the largest sample of size measurements for 
distant early types to date, we have the first opportunity 
to examine the evolution of these galaxies as a function 
of their physical size. To this end we divide our sample 
into three broad size bins as shown in Figure 6, each with 
sample size in excess of 100, spanning the full range of 
typical sizes, namely: (i) 0.5 < r^o < 1.0 h~^'kpc-, (ii) 
1.0 < rso < 2.0 h-^ kpc; and (ill) 2.0 < r^o < 5.0 h-^kpc. 

As we probe to farther distances, the R^p — 24 mag 
redshift reliability cutoff corresponds to ever brighter lu- 
minosities (sec Figure 6). At this magnitude limit the 
completeness of our early-type sample goes to zero. At 
brighter magnitudes the completeness rises rapidly to > 
90% in most redshift intervals. The COMBO-17 survey in- 
cludes a well-defined completeness map constructed from 
Monte-Carlo simulations of the survey and its classifica- 
tion scheme (Wolf et al. 2001, 2003, 2004). Specifically, 
for red-sequence galaxies we have completeness factors as a 
function of iiao for each redshift interval. We have demon- 



strated in §2.4 that the GEMS images of red-selected galax- 
ies are not surface-brightness limited; therefore, the COMBO- 
17 map provides the selection function for our early- type 
sample. For brevity we denote the magnitude and redshift- 
dependent selection function as S{My,z). In each red- 
shift slice of Figure 6 we show schematically the region of 
non-zero completeness (open-ended rectangles). We note 
that the i?ap = 24 mag cutoff cannot be exactly depicted 
in rest-frame F-band owing to the color dependence in 
completeness and the varying transformation between ob- 
served i?-band and rest-frame F-band magnitudes. 

To account for selection effects using distant cluster sam- 
ples, previous studies considered only galaxies brighter 
than a common absolute magnitude cut at all redshifts 
(e.g. Pahre et al. 1996; Schade et al. 1997). However, 
Barger et al. (1998) argued that the consistent analysis of a 
passively-evolving galaxy population must account for the 
evolution of the magnitude limit to avoid underestimating 
the luminosity evolution. 

Our situation is somewhat different in that we wish to 
measure luminosity evolution by direct comparison to the 
well-defined z ~ reference points afforded by the Shen 
et al. (2003) study. Conceptually, we want to passively 
evolve the local sample back to higher redshifts and apply 
the selection function at each redshift. In practice, we de- 
termine the corresponding present-day luminosity for each 
GEMS galaxy of a given size using a Monte-Carlo sampling 
of the mock local distribution that has been subjected to 
the same (but passively faded) selection function as our ob- 
servations. Therefore, for each GEMS galaxy with z, r^o, 
and Mv,z — 5 log^o h, we calculate the luminosity difference 
AMy = Mv,z — ( Afy.o) 1 where {Mv,o) — 5 log^o h is the av- 
erage absolute magnitude from 100 galaxies of the same 
size drawn randomly from the mock z ~ catalog using a 
passively- faded version of S{My, z). We note that this se- 
lection imitates the GEMS galaxy selection while account- 
ing simultaneously for the passive increase in luminosity 
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Fig. 5. — Observed magnitude-size distributions for 928 red-sequence galaxies witli robust c;im2d fits spanning redsliifts 0.2 < z < 1.0 split 
into eight Az = 0.1 bins. In each panel we plot the effective size r.so (geometric mean half-light) and apparent magnitude mgso from Scrsic 
profile fits to the F850LP-band ACS imaging data. We show the expected uncertainties of the fit quantities with error bars at the upper left 
of each panel. We denote n > 2.5 galaxies with solid points and those with n < 2.5 profiles with open squares. The PSF radius (0'.'05) is given 
by the dashed vertical line. Through simulations (see H05) we have established a surface brightness limit of /iso = 24 mag asec^"^ (diagonal 
dot-dashed line) for reliable fits with c;im2d. The additional diagonal lines show surface brightnesses of /150 = 25 (dotted), ^50 = 26 (dashed), 
and = 27 mag asec^'^ (solid), which represent the 85%, 60%, and 0% detection completeness limits, respectively, from simulations of 
71 = 4 profiles (see Rix et al. 2004). We conclude that our early- type sample selection is not surface- brightness limited. For completeness, we 
include the 69 stars misclassified as red-sequence galaxies by COMBO-17 (open stars); we remove these from our further analysis. 
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0.1 1 10 0.1 1 10 

Rest-frame r-band r^Q {h~^ kpc) 

Fig. 6. — Rest-frame luminosity-size relations for 728 early-type galaxies in eight redshift epochs spanning 0.2 < z < 1.0. All galaxies 
have n > 2.5 profiles and red-sequence colors. In each panel we plot the physical size, from gim2d fits to F850LP-band images (corrected to 
rest- frame r-band), against V-hand absolute magnitude. The error bars shown at the bottom right of each panel give the 35% rso uncertainty, 
and the 10% (interpolated for z < 0.7) or 15% (extrapolated for z > 0.7) My — blog^Q h uncertainties. The ridgeline for the z ~ — rso 
relation for early types (n > 2.5) from Shen et al. (2003) is shown in each panel for comparison (diagonal line). The present-day ridgeline is 
the median My — Slog^g h at each rso from our mock catalog (see left panel of Figure 7). The vertical lines delineate the fixed-size bins: (i) 
0.5 < rso < 1.0 /i~^kpc; (ii) 1.0 < rso < 2.0 kpc; and (ill) 2.0 < rso < 5.0 /i~^kpc. The horizontal line attached to the vertical lines 
gives the faint limit of each redshift-dependent completeness map, corresponding to the i?ap = 24 mag COMBO-17 redshift reliability limit at 
each redshift interval. Therefore, the open-ended rectangles outline the non-zero completeness regions of hY,rso, z-space, and the additional 
isolated horizontal line represents the faint limit of the completeness map passively faded to z = 0. We show the physical size corresponding 
to the PSF radius (0'.'05) at each redshift (vertical dashed line). 
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that a present-day galaxy would have at a given epoch. 
Moreover, our method is identical to passively evolving 
(brightening) the z ^ Q mock dataset to the mean redshift 
of each epoch, applying S(My,z), and computing AMy- 

To quantify passive fading, wc adopt the luminosity evo- 
lution of a single-burst stellar population, with formation 
redshift 2;form = 3 and metallicity [Fe/H] = —0.2, using 
the Pegase model (Fioc & Rocca-Volmerange 1997). Be- 
low we discuss other formation redshifts and metallicities, 
including the null no-passivc-cvolution case. In each red- 
shift interval of Figure 6 we show the amount that passive 
evolution shifts the 2; ~ selection function relative to 
S{My, z) by the isolated horizontal line. Since our selec- 
tion function is redshift-dcpcndcnt, it is computationally 
easier to apply a passively-faded version of S{My.z) to 
the z Q mock distribution at each epoch, rather than 
evolving the entire mock catalog at each epoch. 

In Figure 8 we plot the luminosity evolution of the early- 
type population in three separate galaxy size bins. The 
distributions of AAfy values with redshift show evolution 
such that most galaxies of a given size are brighter at larger 
look-back times. The scatter in the data is similar to the 
dispersion given by the log-normal distribution from Shen 
et al. (2003), as shown in the left panel of Figure 7. We fit a 
line (using ordinary least-squares linear regression) to the 
/S.Mv evolution in each panel of Figure 8, constrained to 
intercept the origin {z = 0, AMy = 0). We limit this anal- 
ysis to GEMS galaxies (solid circilcs) within each size bins 
out to a maximum redshift where the selection function is 
> 50% according to the mock 2; ~ catalog for a passively 
fading scenario. Therefore, as outlined with open rectan- 
gles in Figure 6, we do linear fits to galaxy samples from 
the three size bins as follows: (i) 102 with z < 0.7; (11) 227 
with z < 0.9; and (ill) 222 with z < 1.0. We exclude 57 
small (i) and 7 medium (11) size early types (open circles 
in Figure 8). In each panel, the sloped red band gives the 
linear fit to the luminosity evolution with la dispersion. 
The regression dispersion represents the 68-percentile dis- 
tribution of linear fits to 200 bootstrap resamples of the 
(AMy, 2) values. We find significant luminosity evolution 
for early- type galaxies of all sizes, from ^1.0 mag for small 
(0.5 < r5o < 1.0 h~^kpc) systems between z — 0.7 and 
now, to 0.7 mag for the largest galaxies between z = 1.0 
and the present time. We note that the apparently large 
change in Ly — rso slope with redshift found in Figure 6 
is due mostly to our selection function cutting lower lumi- 
nosity galaxies at higher redshift. Nevertheless, once we 
correct for selection effects, the differing amounts of lumi- 
nosity evolution in each size bin of Figure 8 does represent 
a significant change in Ly — rso slope with redshift not 
observed in previous studies. 

In Table 1 we tabulate the slopes of the constrained lin- 
ear fits to the luminosity evolution for each size bin. We 
repeat this analysis using passive fading of S{My, z) that 
is slightly stronger (zform — 2, [Fe/H] = —0.1) and slightly 
weaker (^form = 5, [Fe/H] = -0.5f. We see in Table 1 that 
the measured luminosity evolution is statistically equiva- 
lent (< la differences) between the analyses using three 
different passive fading parameters. In addition, we calcu- 

® In both cases the metallicities are changed in concert with the 
formation redshifts to reproduce the observed color evolution of the 
red sequence; see e.g., Bell et al. (2004b) for more details 



late AMy values under the null hypothesis of no passive 
evolution. As Barger et al. (1998) pointed out, not includ- 
ing passive evolution of the luminosity cutoff (here defined 
by our completeness map) will underestimate the evolu- 
tion. Nonetheless, even under the null condition we still 
find significant luminosity evolution in early-type galax- 
ies of a given size. Furthermore, we find that the amount 
of evolution continues to vary between different fixed-size 
galaxy populations. 

Last, wc compare the constrained fit to the luminosity 
evolution with a fit that is not constrained to go through 
the 2; = 0, no evolution origin. For this test we use the 
same AMy data based on the mock catalog selection with 
-Zform = 3 passive fading assumed, and the same method 
including bootstrap estimation of dispersion. We plot the 
unconstrained fits with light grey bands in each panel of 
Figure 8, and we provide the best-fit linear parameters 
(slope and z = intercept) in Table 2. We point out 
that the intrinsic scatter in the AMy,z correlation, and 
the fairly restricted redshift range over which we fit, both 
combine to produce a less well-constrained fit for the small- 
est galaxies compared with the two larger size bins. In all 
cases, significant fading of the early-type galaxy popula- 
tion at a given size is clearly detected, and the difference 
in fading rates between small and large galaxies is pre- 
served independent of assumptions regarding passive evo- 
lution and fitting method. 

The two methods of fitting luminosity evolution have 
different merits. The constrained fitting is anchored di- 
rectly to 2 = 0, which allows the use of our full knowledge 
of the local universe. Furthermore, the constrained fits are 
quite robust as illustrated by the very narrow dispersion 
we find with bootstrap resampling. On the other hand, 
the unconstrained fits show clearly the differences between 
GEMS and SDSS photometric and size measurement cal- 
ibrations, and for slight deficiencies in our treatment of 
selection effects; the z = offset is < 0.25 mag in all bins. 

3.3. Stellar Mass Size Relation 

We turn our attention now to exploring the stellar mass 
evolution of GEMS early-type galaxies as a function of 
fixed size. In Figure 9 we present the M^ — r^o relations 
for 728 early types in eight Az = 0.1 redshift bins spanning 
0.2 < z < 1.0. As with Ly — r^o, the M* — rso relations 
are quite apparent to at least z ~ 0.8. 

In each panel of Figure 9 we display the M^, — rso re- 
lation ridgeline for z ^ carly-typc galaxies from SDSS. 
Shen et al. (2003) showed that these galaxies also follow 
a log-normal size distribution of the same general form 
as Equation (3) with a stellar mass dependence. These 
authors derived stellar masses from SDSS Petrosian lumi- 
nosities using a model-dependent stellar M^^/L from Kauff- 
mann et al. (2003)^. For our adopted {h is free parameter) 
cosmology, the local M^, — r^o relation has characteristic 
median size 

logio {r5o/h-\pc) = 0.56 logio (M^/iVMq) - 5.52 (6) 
with log-normal dispersion 

0.53 X lO^^M© 



a{M^h^) = 0.29 + 



(7) 



^ Bell et al. (2003) showed that spectrally-derived and color-derived 
stellar masses are the same to within 30%. 
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Fig. 7. — Mock z ~ early-type (ra > 2.5) Lv^.rso (left) and M^,r5o (right) distributions from the Shen et al. (2003) log-normal size 
distributions. In each panel we plot the median (solid line) plus 16-percentile and 84-percentile limits (dashed lines) for the absolute V-band 
magnitude and stellar mass distributions at fixed intervals of r-band half-light physical size. The medians in each panel represent the hy — r^o 
and M* — rso ridgelines that we plot in Figures 6 and 9, respectively. The hy — rso relation has luminosity values that are distributed 
according to the luminosity function of SDSS-selected, red-sequence galaxies from Bell et al. (2004b). Likewise, the — rso relation has 
stellar mass values that follow the early-type, g-band-derived stellar mass function from Bell et al. (2003). 
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Fig. 8. — Luminosity evolution of early-type galaxies at fixed size over the last half of cosmic history. In each panel, we plot the luminosity 
difference AMy between a given GEMS galaxy at redshift z and the average luminosity {{My q) — SlogjQ /i) of a 2 = galaxy with the same 
size, accounting for selection effects (see text), against redshift for all galaxies divided into three size bins as follows: 0.5 < rso < 1-0 kpc (l); 
1.0 < rso 2.0 kpc (11); and 2.0 < rso ^ 5-0 kpc (ill). We fit a line to the data, constrained to intercept the origin (2 = 0, AA/v = 0), 
for galaxies within redshifts limited by the fraction (> 50%) of the mock 2 ~ catalog that we include in our luminosity evolution analysis: 
2 < 0.7 (l); 2 < 0.9 (11); and z < 1.0 (ill). For the linear regression we depict the included (solid) and the excluded (open) galaxies separately. 
The sloped red band represents the best-fit line and its la dispersion from bootstrap resampling of the data (see also Table 1). We repeat the 
linear regression of luminosity evolution for the case where the fit has not been constrained to go through the origin (light grey band) . The 
unconstrained fit results are presented in Table 2. The thin sloped lines are the predicted luminosity evolution of single-age stellar populations 
(computed assuming a Kroupa et al. (1993) IMF using the Pegase code) with different formation redshifts and metallicities: Zf^vui = 2, 
[Fe/H] = —0.1; Zf^,-^ = 3, [Fe/H] = —0.2; and Zform = 5, [Fe/H] = —0.3. The dashed horizontal line represents no luminosity evolution. The 
error bars at lower left of each panel represent the < 15% average uncertainties in COMBO-17 V-band magnitudes. 
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given in observed r-band size (S. Shen 2004, private com- 
munication). Recall that we show in §3.1 that these sizes 
are within 1% of rest-frame r-band. 

Following the method in §3.1, we create a mock cata- 
log of M^-jTso values that represent the M^, — rso distri- 
bution of early-type galaxies at z ~ from Shen et al. 
(2003). Briefly, we construct a 400,000 cell grid of galaxy 
stellar masses and sizes that follow the log-normal dis- 
tribution given by Equations (6 & 7). The grid spans 
8.8 < logio(M*/iVM0) < 12.8 in 800 (0.005 dcx) cells, 
and < rso < 25 h'^kpc in 500 (0.05 /i-^kpc) cells. 
We populate this grid to match the 5-band-derived stel- 
lar mass function for local early-type galaxies from Bell 
et al. (2003). This mock catalog contains over 17.4 mil- 
hon M^,,r5o values. We use this mock catalog to study 
the stellar mass evolution of early types as a function of 
redshift in the next section. In the right panel of Figure 
7, we plot the median value of logj^Q (M*/i^/Mq), with its 
16-percentile and 84-pcrcontilc limits, against r-band size. 
This median provides the Shen et al. (2003) M+ — rso ridge- 
line that we compare to the observed relations at each red- 
shift interval in Figure 9. Here we see that sizes and stellar 
masses of the GEMS early-type galaxies follow a correla- 
tion that is generally consistent with the local M* — r^o 
relation out to ^; ~ 0.8, where the COMBO-17 redshift 
reliability limit begins to cut into the observed relation. 
Under the assumption of simple passive luminosity evolu- 
tion, we expect that galaxies of a given size will maintain 
a constant stellar mass. We investigate this expectation in 
the next section. 

3.4. Stellar Mass Evolution at Fixed Galaxy Sizes 

To directly compare any evolution in the luminosity and 
stellar mass of early types wc must analyze the GEMS 
— rso relation in an analogous manner to the Ly — rso 
relation. As with our luminosity evolution analysis of dis- 
tant galaxies, we calculate how much the M* — r^o relation 
shifts with respect to the local Shen et al. (2003) relation, 
while accounting simultaneously for the selection effects of 
our observations. Just as the COMBO-17 redshift reliabil- 
ity limit at i?ap = 24 imposes an eff'ective absolute mag- 
nitude cut at each epoch, the minimum stellar mass that 
a galaxy has in our sample increases with redshift as ex- 
pected in a magnitude-limited sample (horizontal lines in 
each panel of Figure 9). Therefore, for each redshift bin we 
determine the simple linear relation between M^^ and i?ap 
for each galaxy to convert the COMBO-17 completeness 
map into a rcdshift-dcpendent selection function given as 
a function of logarithmic stellar mass; i.e. ^(logj^Q M^, z), 
where M^ = (M^/i^/M©). We apply S'(logio M^, 2:) to the 
2 ~ mock catalog as before, with one important differ- 
ence - there is no need to correct the selection function for 
the effects of passive evolution because a passively-evolving 
galaxy population is not expected to evolve in stellar mass. 

We divide the sample into the same three fixed-size bins 
and for each GEMS early-type galaxy we calculate the stel- 
lar mass difference Alog^Q M^ — log^Q M^ ^ ~ (logio ^*,o)' 
which is the difference (in log-space) between the stellar 
mass of the galaxy observed at redshift z, and the aver- 
age stellar mass of 100 galaxies of equivalent size drawn at 
random from the z ~ mock M*,r5o distribution, which 
is weighted by ^(logio M^, 0). Therefore, for each GEMS 



galaxy of a fixed size we find the present-day stellar mass 
from the mock catalog using the same selection as our ob- 
servations. The mock catalog is based on a g-band mass 
function; below we discuss the effects that two other mass 
function estimates have on the local sample selection we 
use for this calculation. 

To quantify any redshift evolution of the stellar mass 
of early-type galaxies, we fit lines to the A log;^o ^* 
ues as a function of redshift following the same procedure 
as in §3.2. As before, we limit the analysis to the same 
subset of GEMS galaxies per fixed-size bin (see open rect- 
angles marked i, 11, and iii in Figure 9). In the three 
panels of Figure 10 we plot the linear fits to the stellar 
mass evolution constrained to have no evolution at z = 0. 
The red bands represent the best-fit and the 68-percentile 
distribution of fits (Ict dispersion) from 200 bootstrap re- 
samples. We present the fit results in Table 1. We find 
a modest change in stellar mass for early-type galaxies in 
each size bin. For the two smaller size bins (< 2/i^^kpc) 
the stellar mass appears to have been somewhat larger 
in the past. Conversely, the largest early-type galaxies 
(> 2h~^ kpc) may have had slightly less stellar mass at 
earlier look-back times. We repeat the above analysis us- 
ing local mock Mv,,r5o distributions that follow two ad- 
ditional stellar mass functions based on /C-band (r-band 
concentration-selected) and (/-band {g — r color-selected; 
i.e., Qcoi) from Bell et al. (2003). As shown in Table 1, we 
find that the stellar mass evolution results are independent 
of the mass function we use to construct the present-day 
mock catalog. 

Last, we repeat the linear regression analysis for the red- 
shift evolution of Alogj^QM^ ((/-band-based), but in this 
case we allow the fit to be imconstrained. We present these 
results in the panels of Figure 10 using light grey bands 
(see also Table 2). We see that the unconstrained fits have 
slopes that arc within ^ 2a of fiat; i.e. consistent with 
no stellar mass evolution for galaxies of given fixed sizes 
as expected if they contain passively-evolving stellar pop- 
ulations as suggested by the luminosity evolution (§3.2). 
Furthermore, the intercepts for < 2h~^ kpc galaxies are 
statistically equivalent (i.e., within 3a) to no evolution at 
z = 0, thus, this shows that our stellar mass estimates 
are remarkably well-calibrated with the local SDSS stellar 
masses estimated in an independent manner. We note that 
larger galaxies are about 0.2 dex more massive at 2; = 
than the expectations from SDSS, still a fairly small offset 
given the uncertainties of our estimates. 

4. DISCUSSION 

4.1. Comparison with Previous Work 

Overall the evolution that we see for the Ly — rso rela- 
tion of different galaxy subsets is in good agreement with 
earlier, more-restricted studies based on much smaller data 
sets. The large field of the GEMS survey 15 h^Q„ Mpc 
on a side) permits us to come closer to a cosmologically- 
representative "field" sample of early-type galaxies. Ad- 
ditionally, the large sample size allows us for the first 
time to demonstrate that the strength of luminosity evo- 
lution, at a given size, depends inversely on early-type 
galaxy size: AMy/Az ~ —1.6 ± 0.1 for small galaxies 
(r5o < l/i~^kpc), AMv/Az ~ —1.3±0.1 for medium-sized 
galaxies (1 < rso < 2/i~^ kpc), and AMy/Az -0.7±0.1 
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Fig. 9. — Stellar mass-size relations for 728 carly-typc galaxies with 0.2 < z < 1.0 split into eight redshift slices as in Figure 6. In each 
panel we plot the physical size rso (z-band corrected to r-band) against the COMBO-17 SED-based stellar mass. The error bars shown at 
the bottom right of each panel give the 35% rr,o and the factor of two (0.3 dex) random M* uncertainties. Using our mock M*,r5() catalog 
(see right panel of Figure 7), we show the M* — rso relation ridgeline for z ^ early types (n > 2.5) from Shen ct al. (2003) in each redshift 
slice for comparison (diagonal line). Following the format of Figure 6, we delineate the three fixed-size bins (l, III, III) with vertical lines, 
and we depict the low-mass cutoff (horizontal lines) corresponding to the limit for reliable COMBO-17 redshifts (iJap = 24 mag). Thus, 
the open-ended rectangles outline the non-zero completeness regions of M* , rso , 2-space. The physical size corresponding to the PSF radius 
(0'.'05) at each redshift is shown by the vertical dashed line. 
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Fig. 10. — Stellar mass evolution of early- type galaxies with fixed physical sizes over last half of cosmic history. In each panel we plot 
the stellar mass difference Alogj^gM^ = logjQ ^ ~ (l°gio o) (where = (M*/i'^/Mq)) against redshift for all galaxies divided into 
three size bins shown. A logj^Q represents a measure of the stellar mass evolution from the z ^ Shen et al. (2003) — rso relation (see 
text for details). As with our luminosity evolution analysis (see Figure 8), we perform a linear regression for the stellar mass evolution, both 
constrained (red bands) and not constrained (light grey bands) to intercept the origin (2 = 0, A log^Q = 0). For these fits we use only 
those galaxies (solid circles) within the following redshifts: z < 0.7 (l); z < 0.9 (11); and z < 1.0 (ni). Furthermore, the set of red and grey 
bands represent the 68-percentile distribution of best-fit lines to 200 bootstrap resamplings of the data. We tabulate the fit results in Tables 
1 and 2. The dashed horizontal line represents no stellar mass evolution. The error bars at lower left of each panel represent the 0.3 dex 
random uncertainties in our M* estimates. 



for large galaxies (2 < r^o < 5h^^ kpc). 

Past work has been limited mostly to cluster environ- 
ments (Barrientos et al. 1996; Pahre et al. 1996; Schade 
et al. 1997; Barger et al. 1998; Ziegler et al. 1999), where 
the increased densities provided reasonable sample sizes 
(N < 200) within the small imaging field of WFPC2. 
These studies all found an increase in luminosity with 
look-back time (at a given size) of about 1 magnitude in 
rest-frame i3-band per unit redshift. Allowing for the B 
versus y-band differences and comparing the evolution at 
comparable sizes (~ 3h~^ kpc), we find in comparison no 
evidence for environmental differences in evolution. 

Schade et al. (1999) and Trujillo & Aguerri (2004) have 
selected small early-type galaxy samples {N < 50) with 
WFPC2 observations that are sampled from "field" en- 
vironments at z < 1 in a way that should be statisti- 
cally representative of the cosmic average, and found sim- 
ilar luminosity evolution for galaxies of a fixed size in the 
field as compared to clusters. For example, Trujillo & 
Aguerri found that small early-type galaxies in the HDFs 
were 1.35 ± 0.1 mag brighter in rest-frame T^-band at z ^ 
0.7 when compared to the present-day relation from Shen 
et al. (2003). 

Taking into account the changes in stellar M/L that are 
implied by the color evolution of early-type galaxies in 
the last 8 Gyr, we find that the largely passive luminos- 
ity evolution of the Ly — r^o relation implies little or no 
evolution of the M^, — r^o relation since z ~ 1. Figure 10 
shows that the average stellar mass in early-type galaxies 
of a given size has remained fairly constant over the in- 
terval < z < 1. Concentrating on the constrained fits 
(see also Table 1), we place an upper limit on the change 
in the M^, — r^o relation of 0.25 dex, in the sense that 
galaxies may have been slightly more massive at a given 
size at z ^ 1 than they are at the present. Our findings 
agree with those of Trujillo et al. (2004b), who found very 
little stellar mass evolution since z 2.5 using a sample 



of 168 galaxies of all morphological types from the HDF- 
South. In addition, we find that early-type galaxies of a 
given fixed size had similar stellar masses at z ~ 1 as they 
do today. This is qualitatively consistent with studies of 
the Kormendy relation (e.g., Barger et al. 1998), absorp- 
tion lines (e.g., Ziegler & Bender 1997; Kelson et al. 2001), 
color-magnitude relation (e.g., Kodama et al. 1999; Bell 
et al. 2004b; Holden et al. 2004), and FP (e.g., Kelson 
et al. 2000; van Dokkum et al. 2001; Treu et al. 2002; 
van de Ven et al. 2003; Gebhardt et al. 2003; van der Wei 
et al. 2004), who all found changes in M/L consistent with 
passive evolution of the stellar populations in early-type 
galaxies since z ~ 1. 

The unconstrained fits to the data in Figure 10 may in- 
dicate marginally-significant evolution. In all cases, the 
observed trends within the redshift ranges probed by the 
data are 0.5 dex or less, in the sense that distant early- 
type galaxies may be denser than their low-redshift coun- 
terparts. This trend appears strongest for the 0.5 < < 
1.0/i~^ kpc sample, where we find stronger evolution of the 
Ly — rso relation (Figure 8, left panel). At this stage, it is 
difficult to reliably assess the significance of this possible 
evolution owing to systematic uncertainties of the stellar 
M/L ratios (which are ~ 0.3 dex). In this context, it is 
interesting to note first results of a systematic study of the 
0.7 < z < 1.1 early-type galaxy FP by van der Wei et al. 
(2004); they found that dynamically-derived M/L ratios 
of six intermediate redshift early-type galaxies showed a 
slightly stronger color-M/L correlation than expected in 
a single burst model. The most straightforward interpre- 
tation of this would be an increasing importance of bursts 
of star formation in low mass (therefore small) early-type 
galaxies; aging bursts bias color-based M/L estimates high 
by up to 0.3 dex (Bell & de Jong 2001). Thus, the appar- 
ently strong evolution of the M^, — rso relation for small 
(low mass) galaxies could be enhanced artificially by the 
effect of bursts of star formation, whose frequency is nat- 
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urally expected to be higher at earher times. 

4.2. Understanding the Evolution of Early- Type Galaxy 
Scaling Relations 

Taken by themselves, the observed evolution of the Ly — 

rso and M^, — rso scaling relations is insufficient to con- 
strain the characteristic evolutionary fates of individual 
early- type galaxies. For example, one can imagine mor- 
phological transformations and changes in star formation 
history caused by galaxy mergers, disk re- growth, and fad- 
ing of previously star- forming disks that will cause galaxies 
to drift in and out of the early- type and/or red-sequence 
galaxy populations (e.g., Baugh et al. 1996; Stcinmetz & 
Navarro 2002). Thus, evolution of the scaling relation of 
the populations should be interpreted in terms of the evo- 
lution of the early-type galaxy populations, rather than 
in terms of the evolution of the individual galaxies them- 
selves. 

Observationally, the results of this paper are consistent 
with passive evolution of the early-type galaxy population 
as predicted by the monolithic collapse scenario (Eggen 
et al. 1962; Larson 1974). This model describes the for- 
mation of present-day spherical systems through the col- 
lapse of a massive gas cloud followed by a brief burst of SF 
early in the history of the universe (zform > 2). For this 
model any changes in the observed properties of early-type 
galaxies over time are due to simple passive fading of the 
coeval stellar populations. 

Yet, a large number of works (Chen et al. 2003; Bell 
et al. 2004b; Drory et al. 2004; Conselice et al. 2004; Cross 
c!t al. 2004) have found that the total stellar mass density 
in early-type galaxies, defined either by color or morphol- 
ogy, has built up by roughly a factor of two in the last 8 
Gyr since z ^ 1. Regardless of the mechanisms driving 
this evolution (e.g., mergers, disk fading, etc), the lack of 
drastic evolution in the Mj, — relation indicates that 
early-type galaxies to first order either move along the 
— rso relation as they evolve, or they appear on it when 
they join the sample of early-type red sequence galaxies. 
The lack of strong evolution in the M^, — r^o relation since 
^; ~ 1 is an important constraint that early-type galaxy 
formation theories will have to satisfy. 

It is interesting to note that both disk fading and galaxy 
mergers may naturally satisfy the observational c;onstraints. 
It is possible that some disk-dominated galaxies with rea- 
sonably massive bulges cease to form stars at intermediate 
redshift (owing perhaps to gas consumption or removal of 
its gas supply). As the massive stars in the disk die, the 
disk fades very quickly, increasing greatly the prominence 
of the bulge. Under the assumption that the disk stel- 
lar mass-size correlation does not evolve with redshift (as 
justified by Trujillo et al. 2004b, Barden et al., in prepa- 
ration), and noting that the local disk and bulge stellar 
mass-size correlations are within 0.2 dex of each other over 
a wide stellar mass range (see Fig. 11 of Shen et al. 2003), it 
is quite possible that disk fading would produce early-type 
galaxies that adhere closely to the redshift-independent 
M* — rso relation. 

Another possible formation mode for early-type galax- 
ies is through major galaxy mergers (e.g., Naab & Burkert 
2003; Khochfar & Burkert 2003), where the remnant has 
suffered violent relaxation and is spheroidal and pressure- 



supported (Toomre & Toomre 1972; Barnes 1992). De- 
tailed studies of close pairs have shown that an impor- 
tant fraction of L* galaxies may merge between z ^ 1 
and the present day (Carlberg et al. 1994; Le Fevre et al. 
2000; Patton et al. 2002; Conselice et al. 2003), making 
this a potentially important formation mode for early-type 
galaxies. Yet, it is unclear whether a merger between two 
gas-rich galaxies will lead to a remnant that will satisfy 
our observed lack of significant evolution in the M^, — r^Q 
correlation (sec c;.g., Barnes 2002. for a systematic study 
of gas-rich galaxy mergers; unfortunately the relationship 
between progenitor and remnant size was not explored in 
this work). 

Mergers between gas-poor progenitors are easier to model, 
and numerous studies have found that merging early-type 
galaxies will produce remnants that adhere reasonably well 
to the FP, but slowly drift away from the M* — rso relation 

of the progenitor population (Navarro 1990; Dantas et al. 
2003; Nipoti et al. 2003; Gonzalez-Garcia & van Albada 
2003; Shen et al. 2003). Dantas et al. (2003), Nipoti et al. 
(2003) and Navarro (1990) showed that logio(^remnant/''originai) ^ 
1.21ogio(M 

remnant /Moriginai); that is, ouc dissipationless 
1:1 merger will lead to a size increase of ^ 0.35 dex. For 
a factor of two increase in stellar mass, the M* — r^o re- 
lation (e.g., Equation 6) shows a ^ 0.15 dex increase in 
size. Thus, galaxies undergoing dissipationless merging 
will gradually move towards radii that are larger than 
galaxies not undergoing merging. As outlined earlier, dis- 
sipationless merging is expected to be a much more impor- 
tant process for high-mass galaxies (Khochfar & Burkert 
2003). Therefore, we can focus on the evolution of the 
— rso relation for galaxies with 2 < rso < 5h~^ kpc 
— here the observations place an upper limit of 0.25 dex 
on the evolution of the M* — rso relation zero point since 
2 ~ 1, and show a scatter of ~ 0.5 dex. Therefore, com- 
bining the model predictions (that one major merger will 
move a remnant 0.2 dex from the — rso relation de- 
fined by the progenitor population) and the observational 
upper limits on M* — rso relation evolution and scatter, we 
conclude that the most massive galaxies have suffered at 
most one major dissipationless merger since z ^ 1 on aver- 
age (from the zero point evolution), and that only a small 
fraction of massive early-type galaxies could have suffered 
several major dissipationless mergers (from the scatter). 

4.3. Strength of Luminosity Evolution Dependence on 

Galaxy Size 

It is worth discussing briefly the size dependence in the 
evolution of the Lv — rso relation apparent in Figure 8. 
Focusing on the largest galaxies with 2 < rso < 5 kpc, 
we find that the intercept of the Ly — rso relation is ~ 
0.8 mag brighter at z ^ 1, compared to the present day. 
In contrast, smaller galaxies evolve more rapidly towards 
the present; extrapolated to 2; ~ 1, we find 1.5-2 mag of 
evolution in the Ly — rso relation. 

This difference in evolution is statistically significant, 
thus leaving two possible classes of interpretation: (i) smaller 
galaxies have significantly younger luminosity- weighted stel- 
lar ages, leading to significantly more rapid luminosity 
evolution; and (ii) there is scale-dependence in the evolu- 
tion of the M* — rso relation related to different formation 
routes for low- and high-mass early-type galaxies. At some 
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level, both effects should contribute; we will conclude here 
that we cannot at this stage differentiate between these 
two possibilities. 

The interpretation that smaller galaxies have younger 
stars compared to the larger systems, is consistent with re- 
cent deep determinations of galaxy cluster OMRs at 0.8 < 
^ < 1 (Kodama et al. 2004; de Lucia & et al. 2004), 
where the faint end of the red sequence is systematically 
suppressed compared to the local universe^. This is also 
consistent with analyses of the stellar populations and 
dynamically-derived M/L ratios of morphologically-selected 
early-type galaxies in the local universe (Kuntschner 2000; 
Trager et al. 2000; Thomas et al. 2003), and at z - 1 
(van der Wei et al. 2004), where low-mass early- type galax- 
ies have younger stellar populations than high-mass early 
types. In this interpretation, one expects that the bulk of 
the slope evolution in the Ly — r^o relation is driven by the 
rapid evolution in M*/L of low-mass galaxies, compared 
to their older high-mass counterparts^. 

There is also evidence from studies of local early-type 
galaxies that low-mass early-types have different proper- 
ties than high-mass early-type galaxies (Kormendy & Ben- 
der 1996; Gebhardt ct al. 1996; Faber et al. 1997; Ravin- 
dranath et al. 2001). Stereotypical high-mass early-type 
galaxies have boxy isophotes, steep outer light profiles, 
constant surface brightness (or "cuspy") cores, and are 
supported primarily by random motions of stars in a tri- 
axial potential. In contrast, lower-mass early-type galaxies 
tend towards diskier isophotes, steep power-law light pro- 
files without a resolved core, and derive partial support 
from organized rotation. These differences in properties 
are naturally interpreted in terms of different formation 
mechanisms, where higher-mass early- type galaxies suffer 
from late dissipationless major mergers and lower-mass 
early-type galaxies result from lower mass-ratio interac- 
tions and/or mergers of galaxies with significant gas con- 
tents (Naab & Burkert 2003; Khochfar & Burkert 2003). 
It is also possible that some low-mass early-type galaxies 
could form through disk instabilities in small, very high 
surface density disks (see e.g.. Cole et al. 2000, for one 
implementation of this process). Interpreted in this way, 
the evolution of the slope of the Ly — rso relation could be 
attributed primarily to genuine evolution in the — r^o 
relation, and differences in the rate of M*/L evolution as 
a function of stellar mass may play a secondary role. 

Viewed in this context, the slight differences in stellar 
mass evolution (Figure 10) as a function of size suggests 
that it is impossible to rule out scenario (ii) at this stage. 
Yet, the earlier discussion made it clear that within the un- 
certainties inherent to the estimation of M^,/L from optical 
SEDs, it is impossible to convincingly argue against or in 

^ It is worth noting that one may well expect a change in CMR slope 
with redshift if low-mass galaxies are on average younger. 

^ It is interesting to speculate that owing to biased galaxy formation 
(in the sense that the central parts of high-mass halos collapse at 
higher redshift than the central parts of lower- mass halos), one may 
expect a mild dependence of age on mass, under the assumption that 
one can suppress late cooling in early-type galaxy halos. One may 
expect to see an age difference between cluster and field early-types 
in this interpretation; the observational evidence in this regard is 
controversial, and the most fair statement that can be made is that 
if there is an age diflference as a function of environment, then this 
dependence is rather weak (e.g., Bernardi et al. 1998; Hogg et al. 
2003). 



favor of changes in the — rso relation. Therefore, we 
conclude that disentangling the relative importance of the 
above two mechanisms in driving the redshift-dependent 
slope of the Ly — relation is impossible at this time. 
Future works, utilizing larger samples with spectroscopic 
redshifts and velocity dispersions, will clarify this issue 
considerably. 

5. CONCLUSIONS 

In order to place constraints on the processes driving 
early- type galaxy evolution at z < 1, we have analyzed 
the evolution of the luminosity-size and stellar mass-size 
relations using a sample of 728 early-type galaxies with 
0.2 < z < 1.0 from the GEMS survey in the extended 
CDFS. The sample was selected to have concentrated light 
profiles with Sersic n > 2.5 and rest-frame optical colors 
within 0.25 mag of the red-sequence ridge at its epoch. At 
a given half-light radius r^o, early- type galaxies were more 
luminous in the past. Selection effects were carefully ac- 
counted for throughout using detailed completeness maps 
derived for the parent COMBO-17 photometric redshift 
catalog. We find that the Ly — rso relation has evolved 
in a manner that is consistent with the passive aging of 
ancient stellar populations. In addition, we find evidence 
for an evolving tilt of the Ly — rso relation, in the sense 
that smaller galaxies fade more rapidly towards the present 
day than larger galaxies. Using stellar mass-to-light ratios 
derived from the COMBO-17 SEDs, we rule out any sub- 
stantive evolution in the M* — rso relation, also consistent 
with the simple passive evolution scenario. 

Clearly, these results would be consistent with passive 
evolution of the early-type galaxy population (i.e., passive 
aging of existing stars, no merging), where a younger age 
for smaller galaxies is indicated by the tilting Ly — rso re- 
lation. Yet, this interpretation is too simplistic; a number 
of surveys have demonstrated that the total stellar mass 
in the early-type galaxy population (as defined here) has 
increased by roughly a factor of two since 2 ~ 1. Bearing 
in mind this evolution, our results imply that newly-added 
early-type galaxies follow Ly — rso and M^, — rso relations 
similar to the more established early types. A disk- fading 
origin for the newly-added early-type galaxies appears to 
be consistent with the data. Through comparison with 
models, the role of dissipationless merging is limited to 
< 1 major merger on average since 2: = 1 for the most mas- 
sive galaxies. The predicted evolution for gas-rich mergers 
is not yet robustly predicted, so it is impossible to com- 
ment on this possible formation route. In this context, 
the evolving tilt of the Ly — rso correlation could refiect 
a different origin of low- mass early- type galaxies and/or 
younger stellar populations; the present data is insufficient 
to discriminate between these possibilities. 
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APPENDIX 

Sersic fits 

For our primary analysis we use GIm2d (Galaxy IMage 
2D, Simard et al. 1999, 2002) to fit a single Sersic model 
to the two-dimensional radial light profile of each galaxy 
from the F850LP imaging. The GIM2d software takes into 
account the point-spread function (PSF) by convolving it 
with the best-fit model profile during fitting, and GIM2d 
uses the SExtractor-produced segmentation mask to de- 
blend galaxies from nearby companions. The profile fits 
provide information that allow us to select morphologically 
carly-typc galaxies from the red sequence (as described in 
§2.3.2). In addition, effective size (i.e. half-light radii r^o', 
see §2.4) measurements and other structural parameters 
arc provided by the Sersic fitting. 

Galaxy profile fitting requires a well-defined model of 
the PSF and precise knowledge of the background sky 
level. We use a universal, high S/N PSF derived from 
548 bright but unsaturated stars from the Mark I suite of 
GEMS ACS frames in the F850LP passband (for a first 
impression see Jahnke et al. 2004, plus Jahnkc ct al., in 
preparation). Through detailed testing we find that this 
PSF provides sufficient accuracy for the galaxy fitting (see 
H05). Moreover, with our simulations we explore three dif- 
ferent methods for determining sky values: (1) using the 
SExtractor sky value "local" to each galaxy; (2) adopting 
a fixed constant value per individual ACS frame given by 
the SExtractor "global" estimate; and (3) letting GIM2d 
calculate the mean sky. Wc find that the Sersic fitting 
results are most robust when using the first method, bear- 
ing in mind that other local sky estimates may perform 
equally well or better. 

The Sersic model is represented by seven free parame- 
ters: total intensity, semi-major axis scale length, elliptic- 
ity e, position angle, index n, and model x, y centering. 
For each galaxy fit, gim2d automatically determines the 
initial values and limits for the parameter space to be ex- 
plored by using the segmentation mask, which is affected 
by the choice of detection configuration (see § 2.1). With 
our simulations we find that this method provides reliable 
fits for galaxies with surface brightnesses of fx^o < 24 mag 
asec~^. As can be seen in Figure 5, all red-sequence galax- 
ies arc brighter than this limit. 

We check the fits using a variety of tests. First, we check 
the recovered ellipticity of each fit and find no trends with 
redshift or total mgso magnitude. Likewise, the Sersic n 
distributions show no trends with z or mgso for galax- 
ies with n > 2.5. As expected, we see a slight trend 
such that highly inclined sources have low n values (i.e. 
disk-dominated), and large n values are more often less 



inclined. Lastly, we compare the GIm2d total mgso mag- 
nitudes with the model-independent determinations from 
SExtractor (mag_best). We find a median offset of 0.29 
mag (independent of brightness) with a dispersion of 0.18 
mag, in the sense that GIm2d recovers total galaxy fluxes 
(by integrating the model profile to infinity) that are sys- 
tematically brighter than the SExtractor measurements. 
These tests show that our fits, and the parameters that 
we derive from the fits, arc well-defined. 

For each galaxy, GIm2d finds the best-fit (in a xt sense) 
Sersic model, produces model and residual images, and 
converts each galaxy image into a set of structural param- 
eters including their internal measurement uncertainties. 
In Figures 1 - 3 wc show GEMS F850LP-band galaxy and 
residual postage stamp images for representative early- 
type examples selected randomly from the three fixed-size 
bins out to ~ 1. We note that the majority of galax- 
ies have E/SO visual morphologies, and that an inspection 
of the fit residuals shows that the the fitting does a fair 
job recovering the fiux and overall profile shape. The fit- 
ting results presented here {r^o sizes and Sersic indices) 
will be published in their entirety by Haui31er et al. (in 
preparation). 
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Table 1 

Constrained Linear Fits to Ly — r^o and M^, — rso Relations Evolution 



Ordinate Slopes (m) Notes 







II 


III 




(1) 


(2) 


(3) 


(4) 


(5) 


AMv 


-1.22 ±0.12 


-0.96 ± 0.06 


-0.65 ± 0.09 


unfaded 




-1.65 ±0.13 


-1.33 + 0.06 


-0.72 + 0.08 


Zform = 2, [Fc/H] = -0.1 




-1.63 ±0.12 


-1.28 + 0.06 


-0.71 + 0.08 


Zform = 3, [Fe/H] = -0.2 




-1.59 ±0.13 


-1.26 ±0.06 


-0.71 + 0.08 


zform = 5, [Fe/H] = -0.3 


Alogio (M^/iVMo) 


+0.24 ± 0.05 


+0.12 + 0.03 


-0.13 + 0.03 


g-band-based MF 




+0.27 ±0.05 


+0.14 + 0.03 


-0.10 + 0.03 


gcoi-based MF 




+0.23 ± 0.05 


+0.10 + 0.02 


-0.15 + 0.03 


-band-based MF 



Note. — Wc parameterize the amount of luminosity or stellar mass evolution with the least- 
squares slope from a best-fit linear relation to the given ordinate (1) as a function of redshift, 
constrained to no evolution sd, z — 0. For each galaxy size bin (l, ll, and ill), we give the best-fit 
slopes (2-4), including their la dispersions from bootstrap resampling. The number of galaxy 
measurements used in each fit are A'l = 102, A^n = 227, and Nui = 222. In (5) we list the different 
evolution analyses given by various fadings of the selection function (for luminosity evolution) , and 
by vaxious forms of the adopted stellar mass function (for stellar mass evolution). 



Table 2 

Unconstrained Linear Fits to Ly — rso and M* — rso Relations Evolution 



Ordinate 



Slopes (m) 



Intercepts (6) 



Notes 



(1) 



(2) 



II 

(3) 



in 

(4) 



(5) 



II 

(6) 



III 

(7) 



(8) 



AMv 



Alogio (M^/iVMq) 



-1.60 + 0.46 
-2.56 + 0.56 
-2.40 + 0.54 
-2.38 + 0.52 
+0.69 + 0.20 
+0.63 ± 0.22 
+0.69 + 0.22 



-0.68 + 0.26 
-1.65 + 0.25 
-1.47 + 0.28 
-1.39 + 0.28 
+0.35 ±0.11 
+0.35 + 0.13 
+0.37 + 0.13 



-0.97 + 0.27 
-1.16 + 0.34 
-1.12 + 0.30 
-1.15 + 0.32 
+0.16 + 0.12 
+0.19 + 0.12 
+0.17 + 0.12 



+0.21 + 0.21 
+0.50 ± 0.28 
+0.42 ± 0.28 
+0.43 ± 0.26 
-0.24 + 0.11 
-0.20 + 0.14 
-0.25 + 0.12 



-0.18 + 0.22 
+0.20 + 0.11 
+0.13 + 0.21 
+0.08 + 0.18 
-0.15 + 0.06 
-0.14 + 0.07 
-0.18 + 0.08 



+0.21 + 0.16 
+0.29 + 0.25 
+0.27 + 0.16 
+0.29 ± 0.16 

-0.20 + 0.06 
-0.19 + 0.07 
-0.21 + 0.05 



unfaded 

^form = 2, [Fe/H] = 
2form = 3, [Fe/H] = 
2form = 5, [Fe/H] = 
3-band-based MF 
gcoi-based MF 
if -band-based MF 



Note. — The unconstrained linear fit (least-squares) parameters for luminosity and stellar mass evolution. For each ordinate (1) as a func 

redshift, we give the best-fit slopes (2-4) and intercepts (5-7), including their la dispersions from bootstrap resampling. As in Table 1, we g 
parameters for each galaxy size bin (i, ii, and iii) and each evolution analyses (5). For the unconstrained case we used the same number of 
measurements as with the constrained case. 



